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Systems of tightly knotted, linked, or braided flux tubes will have a universal mass-energy spec- 
trum if the flux is quantized. We focus on a model of glueballs as knotted QCD flux tubes. 
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Plasma physics informs us that linked magnetic flux 
tubes are much more stable than an unknotted single 
loop Q. Linked and knotted flux tubes carry topological 
quantum numbers, and one can think of a knot as a self- 
linked loop. Similar comments apply to braids. 

Our interest will be in tubes carrying quantized flux in 
tight knot and link configurations. If the loops (tubes) 
have fixed uniform thickness and circular cross-section, 
then each knot and link has a completely specified length 
if the configuration is tight, i.e., is of the shortest length 
with the tubes non-overlapping and their cross-sections 
undistorted. If tubes have uniform cross sections, as can 
be approximately the case for many physical systems, 
then the length of the tight knot is proportional to the 
mass (or energy) of the knot. This, we claim, generates 
a universal mass (energy) spectrum for knotted/linked 
configurations of objects of this type. The lengths of 
tight knots were not studied until the mid-1990s |2|, and 
only recently have accurate calculations of large numbers 
of tight knots and links Q become available. These 
results now make it possible to examine physical systems 
and compare them with the knot spectrum. We have 
examined the glueball spectrum of QCD 0,0. Glue- 
balls @ are likely to be solitonic states (See Ref. p| for 
detailed references.) that are solutions to the QCD field 
equations. While QCD will be our main focus, there are 
many more systems where tight knots may play a role. 

In order to decide if a system of flux tubes falls into the 
universal class of having a tight knot energy spectrum, we 
must first investigate the time scales involved. These are 
the lifetime of the soliton t s and the relaxation time r r 
necessary to reach the ground state of a knotted configu- 
ration (i.e., the tight knot state). The soliton lifetime (or 
the corresponding decay width T s = l/r s ) can depend 
on several factors. These include the effects of flux tube 
breaking, rearrangement, and reconnection. The partial 
width for flux tube breaking is non-zero if the produc- 
tion of particle/anti-particle pairs at the break point is 
energetically allowed, for example quark/anti-quark (qq) 
pairs for color electric flux tubes. The partial widths 
can vary widely depending on the particle masses. Re- 
arrangement is a quantum effect where, for example, in 
a linked double do nut arrangement, the loops can tun- 



nel free of each other. Finally, reconnection is an effect 
where tubes break and re-attach in a different configu- 
ration. Such behavior has been seen in plasma physics 
and is of major importance in understanding a variety of 
astrophysical systems. All these processes change topo- 
logical charge, and their partial widths compete more or 
less favorably with each other depending on the param- 
eters that describe the system. 

While no knot lengths have been calculated exactly, it 
is possible to calculate the exact lengths of an infinite 
number of links and braids For links, these calcula- 
tions are possible in the case where individual elements 
of the link lie in planes. For braids, exact calculations 
are possible when the elements of the braid are either 
straight sections or where their centerlines follow heli- 
cal paths. The shortest of all links, the double donut, 
is exactly calculable. The two elements lie in perpen- 
dicular planes and are tori of equal length. The shortest 
non-trivial braid is a helically twisted pair. "Weyl's tube 
formula" |g, Q states that for a tube of constant cross- 
section a normal to a path of length I, the volume of the 
tube in flat 3D is just Vt = Icr. If we have an analytic 
form for the path and a circular cross-section we can find 
Vt ■ This leads to the class of exactly calculable links and 
braids, but since there are no known analytic forms for 
the path of tight knots, their volumes can only be cal- 
culated numerically. We can then calculate or estimate 
the volume and therefore the energy for a corresponding 
physical system. As with knots, the volumes of topolog- 
ically non-trivial tight braids (those where the elements 
are woven together) can only be found approximately. 
While the simple helically twisted braid has a volume 
that depends on the pitch angle which can potentially 
be adjusted by experimental conditions, tight knots and 
links have no such adjustable parameter. 

Let us begin with a discussion of tight links of flux in 
electromagnetic plasma. Movement of fluids can exhibit 
topological properties. Interrelation between hydro- and 
magnetic dynamics may cause magnetic fields, to exhibit 
topological properties as well. For example, for a per- 
fectly conducting fluid, the (abelian) magnetic helicit 
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H = J d 3 x e y '* AidjAk is an invariant of the motion 
and this quantity can be interpreted in terms of knot- 
tedness of magnetic flux lines (The helicity for two 
linked flux tubes with fluxes $i and <& 2 is H = 2n$i$ 2 , 
where n is the Gauss linking number of the two tubes. It 
is straightforward to generalize this to the case of linked 
and/or self-linked thick flux tubes.) A perfectly conduct- 
ing, non-viscous, incompressible fluid relaxes to a state 
of magnetic equilibrium without a change in topology . 
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The system approaches a state of magnetic equilibrium 
by decreasing its magnetic energy by contraction of the 
magnetic field lines. In the case of trivial topology, closed 
curves contract to a point without crossing each other. 
The relaxation eventually leads to a state with zero fields 
(vacuum). If, however, the topology of the initial mag- 
netic fields is non-trivial, the relaxation stops when flux 
tubes are tightly knotted or linked. This happens be- 
cause the "freeze-in" condition forces topological restric- 
tions on possible changes in field configurations and so 
any initial knots and links of field lines remain topologi- 
cally unchanged during relaxation. The energy of a final 
(equilibrium) state is determined by topology. 

By analogy with the abelian case, for a conserved non- 
abelian helicity [l^, we choose the corresponding ex- 
pression for helicity with topological properties, Ti = 
J v (AdA + |^4 3 )- In a perfectly conductive relativistic 
non-abelian plasma, the electric field vanishes in the lo- 
cal frame moving with the plasma. Details are analogous 
to the abelian case. With these facts in mind, we are now 
in position to ask if one could hope to find knotted/linked 
flux tubes in a physical system. 

For several reasons we believe the ideal physical sys- 
tem in which to discover and study tight knots and links 
is Quantum Chromodynamics (QCD). These include: (1) 
QCD is a solidly based part of the standard model of par- 
ticle physics, and much about color confinement and the 
quark model is already well understood in this context, 
making much previous work transferable to the problem 
of tightly knotted flux tubes in this theory. (2) Unlike 
plasmas, fluids or condensed matter systems where flux 
tubes are excitations of some media with many param- 
eters that could hide universal behavior, flux knots in 
QCD can exist in the vacuum. Thus continuum states are 
absent and there are no media parameters to vary and ob- 
scure the universality. Hence, the results in QCD can be 
far less ambiguous. (3) The hadronic energy spectrum 
has been measured over a large range of energies (140 
MeV to 10 GeV) and already many hundreds of states 
are known. We expect that among these, a few dozen can 
be classified as tightly knotted/linked flux tubes states. 
These states must have no valance quarks (i.e., no flavor 
quantum numbers) in order to be classified as glueballs. 

(4) Knotted solitons in QFT are already known to exist. 

(5) One can efficiently search for new glueball states at 
accelerators. (Also, data from older experiments still ex- 
ist and can be reanalyzed to check the predictions of new 
states described below.) 

Consider a hadronic collision that produces some num- 
ber of baryons and mesons plus a gluonic state in the 
form of a closed QCD flux tube (or a set of tubes). From 
an initial state, the fields in the flux tubes quickly re- 
lax to an equilibrium configuration, which is topologi- 
cally equivalent to the initial state. (We assume topolog- 
ical quantum numbers are conserved during this rapid 
process.) The relaxation proceeds through minimization 
of the field energy. Flux conservation and energy min- 
imization force the fields to be homogeneous across the 



tube cross sections. This process occurs via shrinking 
the tube length, and halts to form a "tight" knot or 
link. The radial scale will be set by Aqq D . The en- 
ergy of the final state depends only on the topology of 
the initial state and can be estimated as follows. An ar- 
bitrarily knotted tube of radius a and length I has the 
volume ■na 2 l. Using conservation of flux $>e, the energy 
becomes oc l(tr$ E ) / '(ira 2 ) . Fixing the radius of the tube 
(to be proportional to Aq^ d ), we find that the energy 
is proportional to the length I. The dimensionless ratio 
e(K) = I /(2a) is a topological invariant and the simplest 
definition of the "knot energy" [l^, and can be used to 
fit the correspondence between knot/link energies and 
glueball masses. 

In our model, the chromoelectric fields F^i are confined 
to knotted/linked tubes. After an initial time evolution, 
the system reaches a static equilibrium state which is 
described by the energy density Ee — ^trFoiF 01 — V. 
Similar to the bag model, we have included a constant 
potential energy V needed to keep the tubes at a fixed 
cross-section. The chromoelectric flux is conserved 
and we assume flux tubes carry one flux quantum. To 
account for conservation of the flux, we add the term 
trA{$£;/(7ra 2 ) — n'Foi} to the energy density, where n l 
is the normal vector to a section of the tube and A is 
a Lagrange multiplier. The energy density should be 
constant under variations of the degrees of freedom, the 
gauge potentials A^. This leads to a constant field solu- 
tion, F 0i = ($£/7ra 2 )ni. With this solution, the energy 
is positive and proportional to I and thus the minimum 
of the energy is achieved by shortening I, i.e., tightening 
the knot. 

Lattice calculations, QCD sum rules, electric flux tube 
models, and constituent glue models agree that the light- 
est non-99 states are glueballs with quantum numbers 
J++ = 0++ and 2++ 0. We will model all J++ states 
(i.e., all fj and f'j states listed by the PDG Q), some 
of which will be identified with rotational excitations, 
as knotted/linked chromoelectric QCD flux tubes. We 
proceed to identify knotted and linked QCD flux tubes 
with glueballs, where we include all fj and f'j states. 
The lightest candidate is the /o(600), which we iden- 
tify with the shortest knot/link, i.e., the 2\ link; the 
/o(980) is identified with the next shortest knot, the 3i 
trefoil knot, and so forth. All knot and link energies 
have been calculated for states with energies less then 
1680 MeV. Above 1680 MeV the number of knots and 
links grows rapidly, and few of their energies have been 
calculated (see alternatively, Ref. Q). However, we do 
find knot energies corresponding to all known fj and 
f'j states, and so can make preliminary identifications in 
this region. Our detailed results are collected in Table 
13 where we list fj and f'j masses, our identifications of 
these states with knots and the corresponding knot en- 
ergies. In Fig. ^ we compare the mass spectrum of fj 
states with the identified knot and link energies. Since 
errors for the knot energies in Ref. [2J were not reported, 
we conservatively assumed the error to be 1%. A least 
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TABLE I: Comparison between the glueball mass spectrum and knot energies. 
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/o(600) 


400 - 1200 


2? 


12.6 [4tt] 


768 [766] 


/o(980) 


980 ± 10 


3i 


16.4 


993 


f 2 (1270) 


1275.4 ± 1.2 


2? * 0i 


[6tt + 2] 


[1256] 


/i(1285) 


1281.9 ± 0.6 


4i 


21.2 


1277 






4? 


(21.4) 


(1289) 


/i(1420) 


1426.3 ± 1.1 


5i 


24.2 


1454 


{/ 2 (1430) 


« 1430} d 


5i 


24.2 


1454 + S' 


/o(1370) 


1200 - 1500 


3i * 0i 


(24.7) 


(1484) 


/o(1500) 


1507 ± 5 


5 2 


24.9 


1496 


{/i(1510) 


1518 ± 5} 


5 2 


24.9 


1496 + 8 


/a (1525) 


1525 ± 5 


5 2 


24.9 


1496 + 35 


{/ 2 (1565) 


1546 ± 12} 


r 2 

5i 


(25.9) 


(1555) 


{/ 2 (1640) 


1638 ± 6} 


63 


((27.3)) 


((1638)) 






(2{ * 0i) * 0i e 


[8tt + 3] 


[1686] s 


/o(1710) 


1713 ±6 


6^ 


((28.6)) 


((1714)) 






3i#3l 


28.9 (30.5) 


1732 (1827) 






3i#3i 


29.1 (30.5) 


1744 (1827) 






2J * 2f 


[8tt + 4] 


[1745] 






6 2 


29.2 


1750 






61 


29.3 


1756 






63 


30.5 


1827 






7i 


30.9 


1850 






819 


31.0 


1856 






8 2 o 


32.7 


1957 


/ 2 (2010) 


201ll|° 


7 2 


33.2 


1986 


/ 4 (2050) 


2025 ± 8 


821 


33.9 


2028 






81 


37.0 


2211 






10l61,16 2 


37.6 


2247 


/ 2 (2300) 


2297 ± 28 


818, 9i 


38.3 


2288 


/ 2 (2340) 


2339 ± 60 


9 2 


40.0 


2389 






10i 


44.8 


2672 






Hi 


47.0 


2802 



"Notation n, means a link of I components with n crossings, and 
occurring in the standard table of links (see e.g. 1151 1 on the k th 
place. K#K' stands for the knot product (connected sum) of knots 
K and K' and K * K' is the link of the knots K and K' . 

b Values are from Q except for our exact calculations of 2\ , 2 \ * 0i , 
and (2f*0i)*0i in square brackets, our analytic estimates given in 
parentheses, and our rough estimates given in double parentheses. 

C E(G) is obtained from e(K) using the fit in Figure HI 

^States in braces are not in the Particle Data Group (PDG) sum- 
mary tables. 

e This is the link product that is not 2\ * 2^. 

^Resonances have been seen in this region, but are uncon- 
firmed 0. 

squares fit to the most reliable data (below 1680 MeV) 
gives E(G) = (23.4 ± 46.1) + (59.1 ± 2.l)e(K) [MeV], 
with x 2 = 9.1. The data used in this fit is the first 
seven fj states (filled circles in Fig. in the PDG sum- 



mary tables. Inclusion of the remaining seven (non- 
excitation) states (unfilled circles in Fig. in Table 
where either the glueball or knot energies are less re- 
liable, does not significantly alter the fit and leads to 
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FIG. 1: Relationship between the glueball spectrum E(G) and knot energies e(K). Each point in this figure represents a 
glueball identified with a knot or link. The straight line is our model and is drawn for the fit E(G) = (23.4 ± 46.1) + (59.1 ± 
2.1)e(K) [MeV]. 



E(G) = (26.9 ± 24.9) + (58.9 ± l.O)e(K) [MeV], with 
X 2 = 10.1. Both fits are in good agreement with our 
model, where E(G) is proportional to e(K). Better HEP 
data and the calculation of more knot energies will pro- 
vide further tests of the model and improve the high mass 
identification. 

Knot complexity can be reduced (or increased) by 
unknotting (knotting) operations. In terms of flux 
tubes, these moves are equivalent to reconnection events. 
Hence, a metastable glueball may decay via reconnec- 
tion. Once all topological charge is lost, metastability is 
lost, and the decay proceeds to completion. Two other 
glueball decay processes are: flux tube (string) break- 
ing, which favors large decay widths for configurations 
with long flux tube components; and quantum fluctu- 
ations that unlink flux tubes, which tends to broaden 
states with short flux tube components. As yet we are 
not able to go beyond providing a phenomenological fit 
to these qualitative observations [5j, but hope to be able 
to do so in the future. 



In conclusion, let us return to continuum physics and 
consider a slab of material that can support flux tubes. 
We have in mind a super- fluid or superconductor, but are 
not limited to these possibilities. Assume further that the 
flux tubes carry one and only one unit of flux. Next con- 
sider manipulating these flux tubes. For instance, con- 
sider a hypothetical superconductor where the flux tubes 
are pinned at the bottom of the slab, say by being at- 
tracted to the poles of some magnetic material, and at 
the top of the slab they are each associated with the 
pole of a movable permanent magnet, perhaps a mag- 
netic whisker, or fine solenoid. Assuming the tubes have 
time to relax to tight configurations, the energy released 
should correspond to the universal energy spectrum de- 
scribed above. Another collection of physical systems of 
potential interest are the atomic Bose-Einstein conden- 
sates. For example, laser stirring of dilute 87 Rb atoms 
at 80 nK has produced vortices [lj], which could lead to 
knots and links. 
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